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I. Introduction 


Urethral strictures are a common source of lower urinary tract 
disfunction (Campbell, 1970). Such strictures can occur as congenital 
defects or result from infection, inflammation, or glandular hyper¬ 
trophy. Conventional methods for the diagnosis of urethral strictures 
are time consuming, uncomfortable for the patient and, as they are 
invasive techniques, carry certain risks. 

In the last three years the Academic Hospital in Leiden has 
been evaluating a noninvasive method, called "Urinary Drop Spectro¬ 
metry" (U.D.S.), as an alternative to the conventional techniques for 
initial screening and diagnosis of lower urinary tract obstruction. 

The method of drop spectrometry is based on the observation of the 
external urinary stream by means of an electro-optical instrument. 

Since the external stream breaks up into drops shortly after leaving 
the urethral meatus, the resulting sequence of drops can be transduced 
to a sequence of electrical pulses, in which the statistical properties 
of drop frequency, size, spacing, oscillation, etc. are imbedded. By 
empirical correlation between the statistical properties of the drop 
sequence and conventional diagnostic techniques, the method has been de¬ 
veloped to the stage where gross diagnosis is possible and the effect 
of patient treatment can be traced. For further improvement in the 
diagnostic capability, however, a more fundamental understanding of 
the physiological and hydrodynamic processes is necessary. 


I. A. Purpose of this investigation 


It has been assumed that the variation in the observed U.D.S. 
patterns is due, in part, to an alteration of the turbulent structure 
in the flow downstream from a stricture and thus to an alteration of 
the turbulence in the emerging jet. This alteration in turbulence may, 
in turn, influence the breakup pattern. To test the plausibility of 
this assumption, an investigation has been carried out in two stages: 
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1. R.C.D. Lissenburg has studied the influence of an obstruction 
on the structure of turbulent pipe flow in air. The size and 
shape of the obstruction, the distance downstream from the ob¬ 
struction and the physical characteristics of the canal were 
variables which were investigated. His results will be dis¬ 
cussed in section B of this chapter. 

2. In this investigation, we have scaled the experiments of Lis¬ 
senburg, geometrically and dynamically, to turbulent pipe flow 
of water and have studied the effect of the known alterations 
in the turbulent structure on the process of jet breakup. 

In the present experiments the characteristic dimensions close¬ 
ly approximate the characteristic dimensions of the human male 
urethra; this assures that dynamic similarity to human urination 
is also maintained. The effects of the turbulence have been 
studied by measuring and analyzing the patterns in the drops 
of the disrupted jet. 


I. B. Earlier work 


Despite the many theoretical and experimental investigations on 
the capillary instability of liquid jets, there is not yet a comprehen¬ 
sive theory for the phenomenon. Experimental measurements are most common¬ 
ly presented in the form of a breakup curve (fig. I. 1) which relates 
the length of the continuous portion of the jet to its velocity. The 
breakup curve can be conveniently divided into three regions; a laminar 
region, a transition region and a turbulent region. 

The laminar breakup has been the most intensively studied. At 
low velocities, the jet length is linearly related to the jet velocity 
(laminar region). At higher velocities, the length - velocity relation 
departs from a linear behaviour and, at some critical velocity, reaches 
a maximum (transition region). 
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Rayleigh (1878) laid the analytical foundations for the 
capillary instability of liquid jets through a perturbation analysis 
of the inviscid form of the equations of motion. He showed that only 
axisymmetric disturbances are unstable, and that any axisymmetric 
disturbance whose wavelength was greater than the circumference of 
the jet would grow and eventually lead to breakup. 

The growth rate of the disturbance, and thus the breakup time, was 
found to depend only on the jet dimension and fluid properties and 
was independent of the velocity. As the jet length is equal to the 
product of the breakup time and the velocity, a linear length - veloc¬ 
ity relation is predicted. A critical feature of Rayleighs analysis 
is the dependence of the growth rate on the wavelength of a disturbance. 
A disturbance whose wavelength is equal to 4.51 times the jet diameter 
was found to grow fastest and, providing all disturbances have the same 
initial amplitude (white noise), this most unstable disturbance will 
dominate the jet breakup. The jet should then break into drops sepa¬ 
rated by the distance 


X* = 4.51 D,. 


Continuity of mass 


demands that the drops have the volume and diameter 


v d = i D j 2 * = 3 ‘ 54 D j 3 D d = K89 D j‘ 


The drop frequency can be expressed in terms of the Strouhal number; 


f D, 
Sr - -s-i 

u j 


Since the drop frequency is equal to iL/A*, we have for the Rayleigh 

J 

breakup 


Sr* = D./x* = (4.51) 

J 
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Furthermore, the jet will be unstable to all disturbances with A > tt Dj ; 
i .e. for Sr t < it- 1 . 

Weber (1931) attacked two of the defiencies in Rayleighs analysis. 
He first extended the analysis to include the effect of viscosity, an 
effect treated by Rayleigh only in the limit as the viscosity tends 
toward infinity. Webers result agrees well with experimental measure¬ 
ments in the linear part of the breakup curve and shows that for water 
jets the viscous correction can be ignored. In the second part of his 
analysis, Weber considered the aerodynamic forces on the jet, i.e. the 
effect of the relative motion between the jet and the air. He found 
that the aerodynamic forces enhance the growth rate of axisymmetric 
disturbances and his theory indeed predicts a maximum in the breakup 
curve. The theoretical and experimental values of the maximum length 
and the velocity for which the length is maximum, however, failed to 
agree. Sterling (1975) has modified Webers analysis to account for the 
viscosity of the ambient fluid; the results of his modified analysis 
agree well with experimental values. It was found that aerodynamic 
forces could be ignored provided the Weber number based on the ambient 

density. We = P .4 D j. J_._ , was less than 1. 

a 

Dityakin (1954) has investigated, theoretically, the effect of 
nozzle ellipticity and has shown that if the excentricity (c =vl-b 2 /a 2 , 
a and b are the major and minor axis) is increased, the jet length will be 
shorter and the most unstable wavelength smaller. Thus the drops will 
be smaller and the drop frequency higher. 

Lafrance (1974) has obtained a solution of the equations of motion 
to third order in the perturbations. He has found that the effective 
growth rate of disturbances is very close to that predicted by the 
linear theory and that the drop sizes can be predicted well as a function 
of the wavelength of an initial disturbance. 
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When the jet becomes turbulent, the breakup length is first 
sharply reduced and then increased linearly with velocity. The sudden 
decrease in the jet length is commonly attributed to the large amplitude 
disturbance of the jet surface by turbulent eddies. The study of turbu¬ 
lent jets has been mostly limited to high speed jets ( in application 
to atomization of fuels). There are a few studies on the behaviour at 
low or intermediate Reynolds number and these are largely limited to 
the length - velocity relationships. 

Chen (1964) has found for a turbulent jet, that "the increasing 
size of the surface disturbance as the jet travels downstream is due 
to the propagation of the large scale turbulence distributed within the 
jet. After the influence of the largest scale turbulence is felt, the 
surface disturbance due to the internal turbulence ceases to grow, and 
varicose breakup mechanism continues to desintegrate the jet, producing 
the initial drop. From the orifice to the breakup point, the initial 
disturbance of the jet due to turbulence in the fluid and the surface 
tension forces are thus the predominant factors... 11 . 

Van de Zande (1974) has studied the effect of nozzle length on 
the breakup length of low Reynolds number turbulent jets. With increasing 
nozzle length he has found the transition from laminar to turbulent flow 
is shifting to the region of lower jet velocities. The more fully de¬ 
veloped the turbulence at the orifice, the stronger will be the initial 
disturbance of the jet and thus the shorter will be the breakup length. 

It is clear that, to understand better the effect of turbulence 
on the breakup process, it is important that the turbulent structure 
within the issuing jet is known. This information is provided by the 
study of Lissenburg, i.e. stage 1. 

Lissenburg (1974) has studied the influence of a restriction on 
the structure of turbulent pipe flow of air. He has found that, in the 
turbulent structure of the flow in a tube, the influence of a rigid 
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restriction can be observed at the centreline as far as 40 tube dia¬ 
meters downstream of the restriction and near the wall about 20 tube 
diameters downstream. The shape of the restriction influenced the tur¬ 
bulence such that the more gradually the restriction decreased in the 
downstream direction, the more persistent were the effects. 

The observed effects are most pronounced in the low frequency 
range (Sr < 1). In this range, the turbulent structure does not mono- 
tonically approach the fully developed structure with increasing distance 
downstream from a restriction. Rather, there is first a decrease in the 
turbulent energy until about 10 diameters, then an increase in energy 
to a maximum at about 25 diameters, then again a monotonic decrease to 
the fully developed value. 

Wall roughness obscures the effect of a restriction such that increasing 
the wall roughness of the tube is similar to increasing the distance 
between the restriction and the end cross section. 

If the circular cross section of the tube is brought smoothly to an 
elliptical cross section (designed to simulate the geometric configura¬ 
tion at the downstream end of the urethra, i.e. the fossa navicularis 
and meatus) there is a damping of the turbulent energy in the entire 
wavenumber range. Again the effects are most pronounced in the low 
frequency range. 

He has found further that periodic oscillations induced by vibrat¬ 
ing systems are very persistent and are observed for more than 50 tube 
diameters downstream, even near the wall. Wall roughness, tube curvature 
and the fossa navicularis had little effect on these periodic phenomena. 
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II. Experimental 


The aim of this investigation was to find a possible correlation 
between known alterations in the turbulent structure of a low Reynolds 
number pipeflow and the breakup pattern of the jet. Thus the same test 
conditions have been used here as were used by Lissenburg; i.e. the 
Reynolds number Re^ for the turbulent pipeflow was set to a constant 
value of 5000 and the restrictions used by Lissenburg have been scaled 
geometrically. All experimental variables, except for the modification 
of the test section with the restrictions, were held constant through¬ 
out the investigation. Thus any alteration of the jet breakup must be 
due to the effect of the restrictions. 

The series of experiments carried out can be divided roughly into three 
parts: 


1. The effect on the jet breakup for a given restriction as a 
function of the location of this restriction in the tube, 

2. The effect on the jet breakup of the shape of a restriction, 
placed at a given distance from the orifice of the tube, 

3. The effect on the jet breakup for a periodic disturbance 
caused by an oscillating restriction. 

A flowchart of the experimental procedures is given in fig. II.1. 

For the examination of the jet breakup we have used a small "La¬ 
boratory Drop Spectrometer" (L.D.S.). For a better understanding of the 
phenomena occurring in the jet, high speed cine and a stroboscope were 
also used (section I of fig. II. 1). 

The electrical signal of the L.D.S. has been analyzed either by special 
analog instruments routinely used for patient examination or by a detailed 
analysis on a IBM 370/154 computor (section II). 

The output of the digital computor analysis has been screened and the 
results analyzed statistically (section III). 
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II. A. Apparatus 
II. A. 1. Flowsystem 

The flowsystem consists of a water supply, a flowmeter, a calming 
chamber and a test section 3 a schematic diagram is given in fig. II. 2. 

Mater supply 

The liquid source was an open tank (A) placed three meters above 
the floor. The input of the tank was from the normal water supply. In 
the tank a constant level was maintained by an overflow system (B), two 
meters above the exit of the test section. Because the two meter pres- 
surehead was insufficient, in some extreme cases, a direct connection 
(C) between the output of the tank and the normal water supply was 
provided. 

The flow to the test section was controlled by a valve (D) and measured 
by a variable area, glass tube flowmeter (E). The flowmeter was linear 
over the working range of flow values (fig. II. 3). When use was made 
of the normal water supply by connection (C), it was noticed by obser¬ 
ving the flowmeter that the fluctuation of the flowrate was about the 
same as when the constant headtank was used, (fluctuations were less 
than 0.5% full scale). 

The connection between the calming chamber (F) and the valve (D) was 
made of a flexible tube; this allowed the apparatus to be handled more 
easily. The calming chamber was designed to provide a uniform, laminar 
flow to the entrance tube of the test section. Therefore the diameter 
of the calming chamber was chosen about 15 times the diameter of the 
entrance tube of the test section. Moreover, the upper half of the 
calming chamber was filled with crumpled plastic netting (usually sold 
as pancleaners), while the under half was empty. The two parts are 
divided by a plastic screen of meshsize 1.5 mm (D = 0.25 mm). 
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Test section 

The test section (G) was coupled directly to the calming chamber 
by an entrance tube (I.D, = 5 mm). The test section was held in a ver¬ 
tical position with the jet directed downward. 

When ink was injected into the flexible tube, near the entrance of the 
calming chamber, the water in the under half of the calming chamber 
coloured uniformly. Thus the velocity profile formation in the entrance 
tube of the test section wasn't disturbed by direct streamlines from 
flexible tube to entrance tube. 

Initial observations of the jet formed by this flowsystem showed that 
the turbulence was sometimes intermittent. In order to ensure turbulence 
in the test section a screen of 1.5 mm meshsize (D = 0.25 mm) was placed 
at the beginning of the entrance tube of the test section, while the 
first 25 mm of the entrance tube were roughened by machining V shaped 
grooves in the wal1. 

Precautions were not taken to ensure for fully developed turbulent pipe- 
flow at the most upstream location of the restrictions. Indeed, this 
location is only 23 tube diameters downstream from the calming chamber 
and 10 tube diameters downstream from the 1 mm diameter reduction from 
the entrance tube to the test section. However, as the flowrate was held 
constant for every measurement, for a given location the turbulent ve¬ 
locity profile also was constant and the effect of a restriction on 
this given profile could be measured. 

The test section, calming chamber and most of the water supply were 
made of transparant polyvinylchloride. 

Calibrating the test section 


In order to maintain dynamic similarity between the present 
experiments and the experiments of Lissenburg, the Reynolds number (based 
on the tube diameter) was held constant at Re-^ = 5000. Since 

Gy D 


v 
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and 


0 T = Q/(irD 2 /4) 

we have 

a _ ifDv D „ 

Q — ^ R 0 

Furthermore, in order to maintain dynamic similarity between the present 
experiments and the external urinary stream, we take D = 4 mm. For 
water we have v = 1,10 thus 

0 T = 1.26 m/s 

and Q = 15,7 ml/s. 


This values are within the range of a "normal 11 urination, 
number based on the ambient density, p , becomes 

d 


We 


p a¥j 


0.1 


The Weber 


Thus aerodynamic forces are unimportant. 


It is clear that gravitational forces will affect the breakup process 
through a thinning of the jet as it accelerates; the longer the continu¬ 
ous portion of the jet, the more pronounced will be this effect. In the 
present experiments, however, the continuous jet length was nearly the 
same for all measurements and, as we are looking only for differences 
in the jet breakup with modification of the test section, acceleration 
effects can be ignored. 


II. A. 2. Restrictions and modifications of the test section 


Since the test section was intended to simulate in a very rough 
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way the human male urethra, the L/D ratio of the tube (D = 4 mm) was set 
at 60. The straight open tube simulated the normal male urethra and 
therefore it was the standard configuration to which other configurations 
were compared. Strictures in the urethra were simulated by placing 
restrictions of different size, shape and rigidity at various locations 
in the test section. 

Since Lissenburg has found that a possible effect of a rigid restriction 
is obscured by wall roughness and that this wall roughness has little 
effect on periodic phenomena, no wall roughness has been introduced 
in the test section. 

Rigid restrictions 


In figure II. 4 the five types of rigid restrictions which were 
investigated are shown. 

Restriction 1 simulates a membranous stricture or valve and restriction 
2 a fibrous stricture as does restriction 5. In order to investigate 
whether the effect of restriction 5 was due to the converging and diverging 
part, restriction 3 and 4 were also used. All these restrictions were 
made of a Cu-Al alloyage. 

Elastic restriction and fossa navicularis with elliptical meatus (fig. II. 5) 

An elastic restriction was constructed to test the effect of the 
spontaneous vibrations found by Lissenburg. The shape is almost the 
same as restriction 5. This restriction was constructed by placing a 
specially made rubber tube of 0.05 mm wall thickness and with an inside 
diameter of 4 mm, inside a piece of 5 D length and 4.1 mm diameter 
transparant P.V.C. tube. The connections of the rubber on the plastic 
tube were constructed by turning up the rubber tube ends and glueing 
these to the plastic tube. 

A fossa navicularis with an elliptical meatus, geometrically similar to 
the one used, by Lissenburg, is also tested. This is placed at the end 
of the tube. In the Y-Y plane there is first a widening and then a 
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slight contraction. In the 1-1 plane there is only a smooth contrac¬ 
tion. Together with the elliptical form of the meatus, the fossa 
navicularis can act as a calming chamber and therefore change the 
flow pattern substantially. 

A 5 mm diameter stainless steel tube was used for the test section to 
provide a smooth connection between tube and fossa navicularis. 


II. A. 3. Laboratory Drop Spectrometer 

A schematic review of the L.D.S. transducer and amplifier is 
given in fig. II. 6 and 7. The principle of the L.D.S. is based on 
the interception by a droplet of a light sheet which is directed on 
the sensitive part of a photodiode. 

The almost spherical beam of a small lightsource (A) is converted 
to a nearly homogeneous lightfield by a system of a doublet lens (B) 
and a thin slit (C). A second slit (D) is placed just before the 
photodiode (E) to reduce light reflections from the drops and most 
of the non-parallelism of the lightfield. As the slits have a width 
of 0.2 mm the thickness of the beam is small compared with the diameter 
of a typical drop. Thus the amount of light intercepted by the drop 
(and thus the amplitude of the signal from the photodiode) is proportio¬ 
nal to the local width of a drop in a plane of the light sheet and paral¬ 
lel to the detector (see appendix 1). 

To reduce the pickup of radiated electromagnetic energy, the 
first stage (F) of a four stage amplifier is placed inside the L.D.S. 
and enclosed by a grounded metal box. Also two stages of active filter¬ 
ing at 10 khz reduce both the noise level of sharp spikes from power 
line hash and radiant pickup. 

The negative pulses are about a tenth of a volt superposed on the 7 to 
8 Volt D.C. level of the activated photodiode. The D.C. level is removed, 
the pulses inverted and amplified and the signal is clamped to a zero 
level baseline. A final amplifier provides the adjustment of the output 
voltage. System noise is about 0.5% the height of the standard calibration 
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pulse (see below), which results in a $/N ratio of about 50 dB, in 
practice the noise performance has been limited by the S/N ratio of the 
FM recorder (45 dB). 

Calibrating the L.D.S. 


In calibrating the L.D.S., two aspects must be considered; 

1. The homogeneity and parallelism of the light sheet, 

2. The output of the amplifier related to a constant interrup¬ 
tion of the light sheet. 

re 1) The slit (E) and the photodiode were adjusted against 
each other until the light sheet received by the photo¬ 
diode was homogeneous within 2% over a width of 50 mm. 
This has been measured by dropping from a constant height 
a l inch ball, at several transverse positions within 
the lightfield. The resulting pulses were stored on an 
oscilloscope and compared. If the deviation of the pulse 
height was greater than 2%, the optical components were 
readjusted. 

re 2) The gain of the amplifier was adjusted so that the out¬ 
put pulse from this \ inch ball was 0.50 Volt. The pulse 
height was measured from a stored trace on a calibrated 
oscilloscope. 


II. A. 4. High speed photography 

For a better understanding of the physical phenomena occurring 
as a liquid jet leaves an orifice, two kinds of photography has been 
used; high speed cine with a Hitachi camera and high speed photography 
with a Strobodrum camera. 
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II. B. Summary of experiments 

1. Sensitivity of measurements to experimental conditions. 

A liquid jet and especially a turbulent jet doesn't show break¬ 
up at a unique downstream point. Rather, there is a variation in the 
breakup length. Beyond the initial breakup point, some drops and 
threads of liquid show a secondary drop breakup and the drop pattern 
changes. Thus the minimum distance from the tube orifice for the drop 
pattern to become stable should be known. Once a certain distance for 
recording the drop pattern is chosen, this distance must be kept constant 
for the whole series of experiments. 

To determine this distance, the following experiments have been carried 
out. 


a. The L.D.S. has been located at several distances from the 
tube orifice and recordings of the drop pattern have been 
made. 

b. A high speed cine has been made of the discontinuous part 
of the jet. 

Furthermore the measurements should be reproduceable. Therefore some 
preliminary analyses of drop patterns, recorded several weeks apart, 
have been carried out. 

2. The effect of restriction 5 when placed 10, 20, 30, 40 and 50 

diameters upstream from the tube orifice, 

3. The effect of the different rigid restrictions when placed 10 and 

40 diameters upstream from the tube orifice. 

4. The effect of a self-exciting oscillation of the elastic restric- 

tion. 
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5, The effect of the fossa navicularis and elliptical meatus. 

6. Breakup lengths. 


II. C. General procedure. 


In general, the procedure of the measurements was as follows: 

1) Calibrate L.D.S. system and record four or five calibration 
pulses (0,50 Volt) of the \ inch ball. 

2) Establish a given configuration of the test section. 

3) Position L.D.S. 

4) Adjust liquid flow to get the right Reynolds number. 

5) Record data for 1 minute. 

6) Change configuration of the test section and repeat steps 
2 to 6. 

7) Data analysis. 

re 5) The signal of the L.D.S, system has been recorded on a 

Hewlett Packard 3960 FM Tape Recorder at a recording speed 
of 15 inches per second. Voice annotation was used to label 
the data. 

re 6) To reduce symptomatic errors in the measurements, each ex¬ 
periment was divided into at least 4 parts. For example, 
when the location of the restriction was varied, the se¬ 
quence of measurements was as follows. 

One minute of data segments were recorded, first with no 
restriction in the tube, then with the restriction placed 
10 diameters from the orifice, then 20 diameters from the 
orifice, etc., and finally with the restriction 50 diameters 
from the orifice. The obstruction was then removed and the 
entire sequence repeated, at least 4 times. 
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III. Analysis 


The recorded data, for the most part, have been analyzed digitally 
on an IBM 370 computor. However, sometimes use has been made of special 
analog instruments that are also used for patient examination. These 
instruments analyze only the pulses above a certain constant trigger 
level; the small drops are neglected. The output consist of: 

1) the number of big drops, 

2) the frequency of the drops passing by the lightfield of the 
L.D.S. 

3) a histogram of the time intervals between the big drops. 


III. A, Preliminary digital analysis 

III. A. 1. A.D. conversion 

The signal of the magnetic tape has been converted from analog 
to digital by an IBM 1800 computor. The digitizing rate was set to be 
5,000 Hz and the voltage resolution was about 22,500 units per Volt. 
Before digitizing, the output gain of the FM recorder was checked to 
ensure that the calibration pulses had an amplitude of 0.50 Volts. 

The one minute analog data set for each configuration has been reduced 
to 180 blocks of digital data. Each block consisted of 1024 words and 
at the converting factor of 5,000 Hz, its duration was sec. Each 
set of 180 blocks was preceded by a label block to identify the expe¬ 
rimental run. The digitized data have been written on digital magnetic 
tape. 


Ill. A. 2. Analysis 

The digitized data have been analyzed on an IBM 370 computor. 
The software for these analysis was developed by Prof. J. van Ness as 
part of an earlier investigation of data from the urology department. 
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A flow chart for the program is given in fig. III. 1, while in fig. III. 
2 an outline is given for the treatment of the pulse train from L.D.S. 
to final output. 

In the program, the pulse train is examined and for each pulse above a 
certain trigger level (0.005 Volts), the pulse height, the centroid of 
the pulse, the time since the preceding pulse and the time at which the 
pulse begins and ends has been computed. One data set of 180 blocks 
was divided in 12 groups of 15 blocks, over which groups the means and 
standard deviations for each variable has been computed. Furthermore 
the mean value and standard deviation of the signal has been calculated, 
and the spectra of the beginning and end vector of the group has been 
computed. All these data of one group have been written on magnetic 
tape. For each experiment the mean values and standard deviations 
of the number of drops, the pulse height, the spacing, the signal, the 
spectrum and the histograms of the time interval and pulse height dis¬ 
tribution have also been punched on cards and listed on hard-copy out¬ 
put. 


III. B. Statistical treatment of the output 


Because the number of variables was too large to survey what the 
difference was between two or more configurations, a statistical analysis 
has been carried out consisting of three parts. This should select the 
variables which differentiate best between the configurations which were 
used. 

III. B. 1. Standard mean difference 

Between each pair of configurations the standard mean difference 
for each of the variables has been calculated. 

For n measurements of variable X we have. 

- 1 n 

Mean of X: X = - l X. 

n k=l K 




/, n 

standard deviation S = \ —— e (X, - X) 2 

v n-1 k=l K 


The Standard Mean Difference for X between two configurations i and j 
is defined as 


S.M.D. 



A careful interpretation of all the S.M.D,'s should tell us if there 
is a significant difference in the breakup process for the different 
configuration. A better interpretation can be obtained by taking for 
example the 35 variables with the largest SMD and fed these to some 
statistical discriminant procedures as shown below. 


III. B. 2. Discriminant procedures 

We have used the 35 selected variables and have entered these 
variables for each group of each configuration into a stepwise, discrimi¬ 
nant analysis. This analysis selects, in their order of discrimina¬ 
tion power, the variables which best distinguish between the configura¬ 
tions under consideration. By dividing the total number of groups 
for each configuration into two parts, one for training data, the other 
part for test data, the classification of these test data was our final 
result. If the classification of these test data into the true configu¬ 
ration could be done better by this procedure than by pure chance, 
this was an indication of a difference in breakup pattern. 

For this final analysis we have used two kinds of discriminant 
procedures, a linear and a non-linear. The difference between these two 
methods is the form at the boundary planes between the variables of each 
configuration. For example if there are two configurations the mean 
values for one variable can be calculated for each group. The linear 
procedure has a hyperplane between the two configurations while for the 
non-linear method the curvature of a hypersurfaces can be chosen. 
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For the linear discriminant analysis we have used a program: 
BMD07M Stepwise Discriminant Analysis* Health Sciences Computing Fa¬ 
cility, UCLA. For the non-linear method use has been made of a program 
written by the statists of the Centraal Reken Instituut Leiden. The 
linear discriminant analysis was used to eliminate 25 of the initial 
35 variables and with the remained 10 variables we have used the non¬ 
linear discriminant analysis. For this final analysis all the groups 
of one configuration are training data except one group which have to 
be classified. Doing this several times, each with another group as 
test data, this analysis gives the most important variables and the 
most probable classification of each group for these variables. 
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IV. Experiments and results 


In this chapter the specific procedures and specific results 
will be described for each of the experiments. For all the experiments 
the water temperature was held at 15°C + 1°C. In this temperature 
range the physical properties of water could be regarded as constant. 
Since the water was drawn from the normal water supply the chemical 
properties could also be regarded as constant. The Reynolds number for 
the turbulent pipeflow was set at 5,000 (except for a few special 
experiments discussed later). Any change in the breakup pattern 
from one experiment to another was attributed only to the change in 
the turbulent structure of the pipeflow. 

Survey of experiments 


1) Sensitivity of measurements to experimental conditions. 

2) The effect of restriction 5 when placed 10, 20, 30, 40 and 50 dia¬ 
meters upstream from the tube orifice. 

3) The effect of the different rigid restrictions when placed 10 and 
40 diameters upstream from the tube orifice. 

4) The effect of a self-exciting oscillation of the elastic restriction. 
5} The effect of the fossa navicularis with elliptical meatus. 

6) Breakup lengths. 

IV. A. Sensitivity of measurements to experimental conditions 


In order to know the location where the drop pattern becomes 
stable, the following experiments were carried out. 

IV. A. 1. Measurements 


The distance between the L.D.S. and the tube exit was varied from 
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0.2 m (the breakup point) to 0.42 m in four steps of 0.08 m. At each 
location the drop pattern was recorded for both the open tube at a 
Rep = 7400 and with restriction 1 placed at the orifice of the tube 
at a Re Q = 4200. The data were analyzed by the analog hospital in¬ 
struments . 

In addition to these experiments a high speed cine was made of the 
discontinuous part of the jet from the open tube at Re^ = 7400, 

IV. A. 2. Results 

The data acquired by the L.D.S. and analyzed by the analog 
hospital instruments are shown in fig. IV. 1 and 2, The results of 
the high speed cine are shown in fig. IV. 3. By tracing the path of 
each drop, a history of drops formed at certain place and time has been 
made. This history clearly showed the secondary drop breakup and it 
sometimes appeared that secondary breakup followed drop coalescence. 
However, it was not clear whether this was indeed a coalescence or 
simply an elastic collision between two drops. The history also shows 
a tendency for the drops to travel in clusters. The reason for it could 
be natural or some unobserved pressure fluctuations in the water supply. 
In the further experiments, no such phenomenon was observed either by 
visual inspection of the jet with a stroboscope or by analysis of the 
spectral distribution and the spectrum of the pulse beginnings and 
ends. From the history it can be concluded that from a distance of 
about 0.30 m from the tube orifice the breakup pattern is sufficient 
stable. In the spectral distributions, however, the difference in the 
drop pattern becomes small only for distances greater than 0.6 m from the 
tube orifice; at this distance the breakup pattern is stable. The dif¬ 
ference between fig. IV. 1 and fig. IV. 2 could be expected because 
restriction 1 placed at the tube orifice, acts as an orifice with dia¬ 
meter half that of the open tube and therefore produces smaller drops 
and a faster stream than the open tube does. 

It was determined, mostly from the high speed cine, that a distance of 
about 0.45 m from the tube orifice would be sufficient to provide an 
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almost stable drop pattern. 

The reproduceability of the measurements has also been tested. In some 
preliminary experiments a few recordings of the drop pattern from the 
tube without restriction were made, several weeks apart. There were 
small but monotonic changes in the patterns with time. To avoid the 
complications which such "aging" effects would introduce into the 
analysis, all experiments for which drop patterns were compared were 
completed within two hours. 


IV. B. The effect of restriction 5 when placed 10, 20, 30, 40 and 50 

diameters upstream from the .tube orifice 


For this experiment we have chosen restriction 5 because Lis- 
senburg has found for this restriction still a change in the turbulent 
structure at the orifice- when the restriction was located at a distance 
of 40 diameters upstream from the orifice. 


IV, B. 1. Measurements 

Of each configuration i.e. the restriction located at a given 
distance from the tube orifice, six recordings of the drop pattern 
have been made. The sequence was as follows: open tube, L/D = 10, 20, 
etc., 50; open tube, L/D = 10, etc. 

IV. B. 2. Results 


The data obtained for the drop pattern with restriction 5 placed 
at 10, 20, 30, 40, and 50 diameters upstream from the orifice are given 
in figures IV. 4 and 5 and tables IV. .1 and 2 together with the data ob¬ 
tained for the open tube. 

Clearly, a variation in the location of the restriction does 
not lead to large differences in the drop pattern. Thus the much 
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greater differences in the pattern observed with patients cannot be 
attributed simply to a variation in the location of an obstruction. 

The pulse height distribution (fig. IV. 4) has a typical bimo- 
dal form, i.e. the discontinuous stream consists of large "main" drops 
and small "satellite" drops. Whereas they are observable differences 
in the distributions, there are no monotonic or systematic changes 
with the location of the restriction. The data are, however, nearly 
bounded by the distributions for L/D = 10 and L/D = 30; the drops for 
L/D = 30 are smaller and the mode for the satellite drops is narrower. 
This agrees qualitatively with the results of Lissenburg that the ratio 
of low frequency to high frequency turbulent energy is a minimum for 
L/D = 10 and maximum near L/D = 30. 

The shape of the time interval distribution (fig. IV. 5) 
suggests that the drop arrival times are Poisson distributed. The 
time interval distributions for all configurations are nearly identi¬ 
cal . 


The plots of the spectra of the pulse beginnings and ends 
vector aren't given because the differences for the curves of different 
configurations are small and no visual information can be drawn. The 
same holds true for the time interval distributions. 

The change in every variable with a change in the location 
of the restriction is small (table IV. 1); there is no single variable, 
which clearly distinguishes the location of the restriction. That is, 
the standard mean differences for all variables and for all configu¬ 
rations are small. Thus it is necessary to apply a multivariable analysis 
to see if systematic differences exist. The results for both the linear 
and non-linear discriminant analyses are given in table IV. 2, 
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Linear discriminant method 


Three matrices are given A, B, and C. For A and B the same 
variables were used. Matrix A is the classification after the selection 
of the 10 variables which best distinguish, while for matrix B all 
the 35 variables were used. The difference between the two matrices 
can be explained by the number of variables, 35 variables are too 
many for only 48 groups. 

Within the 10 best variables there are 5 from the pulse height distri¬ 
bution, so an analysis is carried out with variables chosen only from 
the pulse height distribution; matrix C. These three matrices show, 
in the diagonal, a sudden decrease in correct classification when the 
restriction was placed more than 20 diameters upstream. At this and 
greater distances the different patterns are difficult to separate. 

That is, it is difficult to determine if the restriction is at L/D = 

30, 40, or 50 but it is more clear that the restriction is not at L/D - 
10 or 20. 

Non-linear discriminant method (matrix D) 


When the 10 best variables of the linear discriminant analysis 
were entered in a non-linear discriminant method, two variables from the 
spectra of the pulse beginnings and ends vector are first chosen followed 
by a variable from the pulse height distribution. A result, not fully 
explained, is the much better classification of the patterns for the 
location 20 and 50 diameters upstream from the orifice. The matrix A 
shows the same but smaller effect. The matrix D shows that only the 
pattern for the location L/D = 20 is unique, while the patterns for 
the other locations look, for the most part, like the patterns for the 
location L/D = 50. 
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IV. C, The effect of the different rigid restrictions when placed 10 
and 40 diameters upstream from the tube orifice 


The aim of this investigation has been twofold: 

1) the effect of different shape and size of a restriction on 
the breakup pattern. 


2) comparison between the restrictions for the persistancy of 
a possible change in breakup pattern, when the distance 
between the location of the restriction and the tube orifice 
has been increased. 


IV. C. 1. Measurements 

At both locations four recordings of the breakup pattern have 
been made for each restriction. The recording sequence was as follows: 
open tube, restriction 1 placed 10 diameters upstream from the tube 
orifice, restriction 1 placed at 40 diameters upstream, restriction 1 
placed 10 diameters upstream, etc. (four times), then open tube again 
followed by restriction 2 placed 10 diameters upstream, etc., until 
restriction 5. 

For restriction 1 the two-meter pressure head was to low to provide 
a Reynolds number of 5000. Only for this restriction the recordings 
of the breakup pattern have been made while the desired flow rate was 
supplied by the normal water supply. 

The variables that were used in the statistical analysis were for 
both locations the same. They have been gathered from the 
analysis of the data of the configuration where the restrictions were 
placed 10 diameters upstream from the tube orifice. 


IV. C. 2. Results 


The data obtained for the drop patterns of the different res¬ 
trictions are given in fig. IV. 6 and 1 and in table IV. 3, 4 and 5 
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when the restrictions were placed 10 and 40 diameters upstream from 
the orifice, respectively. 

An overall comparison between the two pulse height distributions indi¬ 
cates a slight shift to the larger drops for the L/D = 10 curve compared 
with the L/D = 40 curve of the same restriction. When placed at a 
distance L/D = 10 from the orifice, the restrictions show differences 
in the pulse height distribution, which are greater than observed with 
restriction 5 in different locations. 

Translating the differences in the curves to the effect of different 
size and shape of a restriction on the drop pattern is difficult. In 
general, it can be said that any restriction placed at a distance of 
10 diameters upstream from the tube orifice has a specific influence 
on the breakup pattern. When restriction 2, which has a greater 
throat than the other restrictions is not considered, the influence 
of the upstream or downstream shape can be deducted as follows (fig. 

IV. 6): 

If we assume, the downstream shape of a restriction is important, 
then restriction 1 and 3 should have almost the same effect on the 
drop pattern. The pulse height distributions show that this is not 
true. Thus the upstream shape can be important or a combination of 
upstream and downstream shape. If we assume, the upstream shape is 
important, then restriction 1 and 4 should have almost the same effect. 
The pulse height distribution indicates that this is nearly true. 

The curves of these restrictions are more similar than the curves of 
restriction 1 and 3. Thus the upstream shape of a restriction appears 
to be the dominant factor for the breakup pattern. The pulse height 
distribution for the restrictions when placed at a distance of 40 dia¬ 
meters from the tube orifice (fig. IV. 7} confirms the conclusion made 
above: the curves for the restrictions 3 and 5 are almost the same 
but are different from the curves for the open tube. The curves for 
restrictions 1 and 4, on the other, are nearly the same as for the 
open tube. Thus, the upstream shape leads to the most persistant 
effects. 

In table IV. 3 the mean values andstandard deviations of some basic 
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variables are given. Again, there is no variable that can distinguish 
the different restrictions except restriction 2. In table IV. 4 and 5 
the results of the discriminant analyses are given. Looking at the 
distribution of the misclassified groups, one notices that restric¬ 
tion 1 and 4 look like each other as do restriction 3 and 5. 

The effect of the larger throat of restriction 2 seems to be a 
production of larger drops. 

IV. D. The effect of self-exciting oscillations of the elastic res- 
triction 


Lissenburg found in his investigation of an elastic restriction; 
that under certain circumstances this restriction exhibited self-exciting 
oscillations. Similar sets of experiments has been done in this investi¬ 
gation. 

IV. D. 1. Measurements 

When the external pressure on the membrane of the elastic res¬ 
triction had reached a value of about 6.10 3 N/m 2 , a self-exciting 
oscillation of the membrane occurred. This could easily be observed 
by looking at the signal of the drop pattern on the oscilloscope. In 
addition a faint low frequency tone could be heard. 

However, the whole system of oscillating elastic restriction and flow 
system wasn't stable. The amplitude of the oscillation decreased with 
time until the effect on the jet was no longer observable. After several 
attempts, one recording of three minutes duration of the drop pattern 
was made in which the phenomena were stable; the elastic restriction 
was 40 diameters upstream from the tube orifice. In addition to this 
stabilization problem, there occurred also a construction problem; 
no two elastic restrictions could be constructed to give the same 
properties. 
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IV. D. 2. Results 


The pulse height distribution for the oscillating elastic res¬ 
triction placed 40 diameters upstream from the tube orifice is given 
in fig. IV. 8. In the same figure also a distribution for the open 
tube is plotted. The distribution for the oscillating elastic res¬ 
triction also shows a pattern with two modes. These modes are more 
slender than those of the open tube, thus the jet breaks up more 
into preferential drops. This is especially pronounced in the mode 
for the larger drops. 

The oscillation of the flow field forces the jet to break 
up in smaller "bigger drops" (the shift to the left of the bigger drops 
mode). Thus an oscillation of the flow field leads to such a strong 
initial disturbance at the tube orifice that, even if this gives a 
wave length much smaller than the most unstable one (the Rayleigh 
wave length), this disturbance controls the jet breakup pattern. 

IV. E. The effect of the fossa navicularis and elliptical meatus 

Lissenburg has found that the effect of the fossa navicularis 
and elliptical meatus is to damp the turbulent energy across the entire 
wave number space. 

It has been assumed, therefore, that the fossa navicularis and ellipti¬ 
cal meatus on the jet breakup pattern should show a same effect. Dif¬ 
ferences in drop pattern caused by differences in restrictions or lo¬ 
cation of the restriction should stay the same therefore or should be 
less. Thus only the effect of the fossa navicularis and elliptical 
meatus on the drop pattern of the tube without restrictions has been 
studied. This has been done in combination with an examination of the 
influence of the jet velocity on the drop pattern. 

IV. E. 1. Measurements 


First, a five minutes recording of the drop pattern was taken 
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with the Reynolds number for the 5 mm diameter tube set at 4000. Then 
the fossa navicularis and elliptical meatus were coupled to the tube, 
the L.D.S. was lowered to compensate for the added length of the fossa 
navicularis and elliptical meatus. Again a five minute recording was 
taken of the drop pattern of this configuration; the Reynolds number 
was the same. This sequence was then repeated for Re Q = 5500 and 
Re^ = 8500. 

The drop pattern from the fossa navicularis and the elliptical meatus 
showed a significant difference from those of the open tube. This 
difference could be observed on the oscilloscope. The number of pulses 
increased and the pulses showed a more uniform shape and size. 

Also was it noticed, by visual and stroboscopic inspection, that the 
breakup length of the jet shortened remarkable. 

IV. E. 2. Results 


The pulse height distribution of the drop pattern for the open 
tube cases are given in fig. IV. 9. The three curves all show a bimodal 
pattern, but with increasing Reynolds number the smaller drops mode 
becomes less important. The top of this mode decreases and shifts to 
the left and while the base stays almost the same the mode becomes 
asymmetric. The mode of the bigger drops shows, with increasing Rey¬ 
nolds number, a flattening of the top while the maximum first decreases 
and then increases. The base of this mode first shifts and then 
widens. Thus, with increasing Reynolds number, the preferential size 
of the satellite drops decreases while the number of these small drops 
show a strong decrease. The preferential size of the big drops 
stays the same. However, the preference of the jet to break up into 
drops of the same size decreases (flattening of the mode). 

The pulse height distribution for the fossa navicularis and 
elliptical meatus cases are given in fig. IV. 10. With increasing 
Reynolds number the bimodal pattern almost changes into an unimodal of 
only primarily drops. The two modes are almost symmetric. The top of 
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the smaller drops mode decreases and shifts to the right while the 
base stays the same. For the bigger drops the top of the mode also 
decreases but there is only once a shift of the top (Re^ = 4000 -+■ 
5500). The base increases with increasing Reynolds number. Upon 
comparing the two figures we can see that the effect of the fossa na- 
vicularis and elliptical meatus on the drop pattern consists of: 

1) Producing a smaller size of the preferential drop for both 
the smaller and the bigger drops mode. 

2) Decreasing the number of smaller drops and increasing that 
of the bigger drops. 

3) Transferring the asymmetric modes of the open tube drop 
distributions to almost symmetric ones. 

IV. F. Breakup length 


A rough measurement of the breakup length of the jet was made 
for each of the experimental conditions described above. The distance 
from the tube orifice to the most probable breakup point was measured 

by a visual inspection of the jet with the stroboscope. For all the 

experiments with a flowrate of 15.7 ml/s the breakup length of the jet 
was 0.19 m + 11. This result (a constant breakup length) confirms the 
assumption that the gravity will not have a different effect on the dif¬ 
ferent experiments. For example, if the breakup length for one experi¬ 
ment was much shorter or longer then there should be a gravitational 
influence on the breakup pattern due to the difference in jet diameter 
at the breakup point. 

For the experiments with variable flowrate, the breakup length 

increased with increased Reynolds number. The following breakup lengths 

were measured (I.D. = 5 mm; L/D “ 60). 
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Rep 

breakup length 

4000 

0.21 m 

5500 

0.25 m 

8500 

0.38 m 


These values have an almost linear relation, which indicates that for 
this experimental setup the jet lies in the turbulent region of fig. 

I. 1. 


The effect of the fossa navicularis and elliptical meatus was 
a decrease in the breakup length, compared with the open tube. The 
relation between breakup length and Reynolds number again is linear. 

Rep breakup length 


4000 

5500 

8500 


0.17 m 
0.21 m 
0.32 m 
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V. Succeeding investigations 


When the effect of different restrictions on the jet breakup 
was studied (chapter IV. C), the analysis of the recorded data showed 
a significant difference for two of the five investigated restrictions, 
for both the locations L/D = 10 and 40 diameters distance from the tube 
orifice. Restrictions 2 and 4 showed an anomalous behaviour, especial¬ 
ly in the pulse height distribution. 

As Lissenburg did not investigate the effects of restriction 3 and 4 
on the pipe turbulence, these measurements were taken as a part of 
this study. The apparatus and procedures were identical to those 
used by Lissenburg. The results are given in fig. V. 1 and 2 for 
the mean velocity and turbulent intensity profiles and the spectra 
of the turbulent energy for both the open tube and restriction 3 
and 4. 

In these experiments,no anomalous behaviour was found for restriction 
4 compared with the open tube case. The conclusion of Lissenburg, 
that the more gradually a restriction diverges in downstream direction 
the more persistent the effects, has also been confirmed. 

A further study of the anomalous behaviour of restrictions 2 and 4 
(when water was used) revealed that the anomalous data occurred only 
in isolated segments of the recordings and that the isolated segments 
were sequential. That is, this behaviour was nor occurring intermit¬ 
tently but rather occurred only over one portion in each of the ex¬ 
periments with restriction 2 and 4. It was not possible to reproduce 
this behaviour in later experiments. 

The anomalous data have been analyzed separately and the 
pulse height distributions data are given in fig. V. 3,together with 
the pulse height distribution for the self-exciting oscillation of the 
elastic restriction. From the similarity of the curves it is concluded 
that a spontaneous oscillating flow field is possible downstream from 
a restriction and that this oscillation can effect a change in the 
jet breakup. For the elastic restriction an external pressure of 
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6.10 3 N/m 2 was needed to effect a same behaviour of the jet. It is 
assumed, therefore, that such an anomalous behaviour of the jet is 
due to the existence of a more or less discrete eddy, persisting 
with distance downstream from the restriction. A physical reason for 
the existence of such an eddy can be the following. Assume a two 
dimensional system. When the divergency of a restriction in down¬ 
stream direction is greater than the angle of the "emerging jet 11 
leaving the throat of the restriction, it is possible that the 
"emerging jet" tends to lie on at one side of the restriction. This 
will be caused by the under pressure in the liquid layer between 
"emerging jet" and the wall of the restriction. However, this is not 
a stable flow because for the same reason the "emerging jet" tends to 
lie on at the other side of the restriction. Thus an oscillating flow 
field is possible. In the three dimensional case it is also possible 
that in addition to this oscillation a rotation of the flow field 
occurs, which can lead to a very persistent eddy. 
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VI. Conclusions 


First, the specific conclusions for each set of experiments 
will be drawn. Then the combined results from stage 1 and 2 will be 
discussed and final conclusions will be drawn. 

VI. A. Specific conclusions 

VI. A. 1. The effect of restriction 5 when placed 10, 20, 30, 40 and 
50 diameters upstream from the tube orifice 


A restriction located at a distance to the orifice less than 
25 diameters has an influence on the jet breakup. Beyond this distance 
the effect rapidly decreases. However, at the location L/D = 50 there 
again appears to be a small effect. Because the effect on the drop 
pattern does not decrease monotonically with increasing distance of the 
restriction from the tube orifice, the localization of the restriction 
in the tube will be restricted to a qualification: The restriction 
has been placed at a distance to the tube orifice less than or greater 
than 20 diameters. For this qualification, the pulse height distribu¬ 
tion provides the most important variables. 

VI. A. 2, The effect of the different rigid restrictions when placed 
10 and 40 diameters upstream from the tube orifice 


The differences in breakup pattern for the different restric¬ 
tions depend strongly on the location of the restriction. When the 
restrictions are placed at a distance of 10 diameters upstream from 
the tube orifice each restriction had his own effect, more or less. 

The restrictions with a diverging shape in upstream direction influence 
the breakup pattern most. At a distance of 40 diameters upstream from 
the tube orifice, the patterns for the different restrictions become 
almost the same. But still the two restrictions with a smooth diverg- 
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ing shape in upstream direction have most difference with the open 
tube. 

VI. A. 3. The effect of self-exciting oscillation of the elastic 
restriction 


Evidently, a self-exciting oscillation influences the jet 
breakup pattern. The oscillation causes the jet to break up into small 
drops, while the distance of the oscillating restriction from the ori¬ 
fice seems of little importance. 

From the anomalous behaviour of the restrictions 2 and 4, it is evident 
that the initial oscillation can be very faint and yet control the jet 
breakup. 


VI. A. 4. The effect of the fossa navicularis and elliptical meatus 


The fossa navicularis and elliptical meatus, when attached to 
the open tube, alters the pulse height distribution from a bimodal 
(large "main 11 drops and small "satellite" drops) to a nearly unimodal 
distribution. This trend is more pronounced as the Reynolds number 
is increased. The main drop size is decreased; this is primarily an 
effect of the elliptical meatus. 

VI. A. 5. Analysis of the recorded data 


For investigations carried out on the breakup pattern of a 
liquid jet at a constant flowrate, the drop size distribution carries 
the most information. 

To interpret the differences in breakup pattern for the different res¬ 
trictions or locations a statistical analysis is necessary. (The differ¬ 
ences in percentage correct classification for especial the experiments 
with changing location, can be the result of the differences in the two 
statistical analysis methods, which differences can have great effect 
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when the cases to classify are more or less equal. A further study 
of this problem is necessary.) 


VI. B. Comparison between the experiments with air and water 


The results from the experiments carried out in air show that 
some restrictions can be observed, at the centreline, as far as 40 dia¬ 
meters downstream from that restriction. 

The experiments carried out with water still show some effect on the 
drop pattern when the restriction was placed at a distance of 50 dia¬ 
meters from the tube orifice. These two distances of 40 and 50 diameters 
have been the longest distance from the tube orifice for which location 
measurements have been taken. 

The results for the influence of the size and shape of a restriction, 
however, do not agree between stage 1 and 2. In air the effects are 
more persistent the less the restriction diverges in downstream direc¬ 
tion. For the experiments with water the restrictions which diverge 
in the upstream direction have the most persistent effects. The down¬ 
stream geometry appears to have little effect. 

When there is an oscillating source in the air flow field, 
the turbulent energy is transferred to the higher frequency eddies with 
distance downstream. For the experiments with water the pulse height 
distribution tends to shift the top of the main drops mode towards 
the bigger drops, with distance downstream. 

In the experiments with air the fossa navicularis with elliptical mea¬ 
tus damps out the turbulent energy across the entire wave number space 
but this effect is especially pronounced at the low frequency range. 

For water, the effect on the jet breakup is a more uniform distribution 
for the bigger drops. These bigger drops, however, are smaller than 
those of the tube without fossa navicularis and elliptical meatus; 
this is primarily an effect of the elliptical meatus. 
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When the turbulent energy spectra are compared with the pulse 
height distribution, especially for the experiments carried out on the 
fossa navicularis and the oscillating restriction, the following conclu¬ 
sion can be drawn. When the energy in the lower frequency range de¬ 
creases the small drops become less in number (fossa navicularis) and 
when the energy in the lower frequency range increases, the drops be¬ 
come smaller (oscillating restriction). Moreover when one certain 
eddy in the lower frequency range has sufficient energy, this eddy 
will control the drop pattern (anomalous jet breakup). 

This assumption is supported when the pulse height distribution and 
the energy spectra are compared for the experiments where the distance 
to the orifice has been increased (fig. IV. 4 and fig. V. 2). The 
energy spectrum shows, in the low frequency range, a decrease in 
energy compared with the open tube when the restriction has been placed 
within a distance of 10 diameters from the tube orifice. A maximum in 
energy for the low frequency range was reached at 20 diameters. 

In the pulse height distribution the location L/D = 10 shows for both 
the right flank of the smaller drops mode and the left and right flank 
of the bigger drops mode a small shift to the right, compared with the 
open tube. For the location L/D = 20 only the left flank of the bigger 
drops mode shows a remarkable shift to the left. 

VI. C. Final conclusions 


The jet breakup pattern is influenced by a restriction located 
upstream from the orifice. The influences of the size, shape, and lo¬ 
cation, however, lead to much smaller differences in the drop pattern 
than is observed in the external urinary streams of human males, and can 
not be regarded as the dominant influence. 

There exists a correlation between the results of the measure¬ 
ments with air and water. An alteration in the turbulent energy of the 
eddies, especially in the lower frequency range, causes a specific 
reaction in the breakup pattern. 
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A diverging shape in upstream direction of the restriction 
causes more effect on the jet breakup than does a restriction with 
a sharp edge in upstream direction. 

The effect on the jet breakup caused by a restriction with 
a diverging shape in upstream direction is more persistent with distance 
downstream than the effect caused by a restriction with a sharp edge 
in upstream direction. 

Periodic oscillations change the breakup pattern totally as 
does the fossa navicularis and elliptical meatus. However, the change 
is not the same. 

To observe differences in the jet breakup pattern at a constant 
flowrate, the pulse height distribution is most important. The time 
interval distribution is of only minor interest. An intensive statis¬ 
tical treatment of the data is necessary to detect the effects of size, 
shape, and location of a restriction. 
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APPENDIX 1 


Basic principles of the optical-electro transducer 


The pulse height, caused by an interruption of the light sheet by a 
drop is linear with the amount of light intercepted. When the drops 
are assumed spherical the following relations can be calculated: 

pulse height h (time) 

width of the interruption 1 (time) 

calibration factor cf 



drop velocity U 
pulse width tw 


l 

\ 
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t = - ^ + x = 0 


t=+^^x=d 


and 


tw = 


U 


x = (t + ^) U 


substitute in (2) 


i (t ) = d [i -(^i) 2 p 


substitute (5) in (1) 


(t) = cf.d [l - (^) 2 p 


= 0 


t S 


else 


< d 


7U 


(3) 

(4) 

(5) 

( 6 ) 


Only for spherical drops the drop velocity can be calculated from the 
pulse 


for t = o:h(o) = h m=v = cf.d 
v ' max 

Substitute (3) leads to: 

h 

I! _ max 

~ cf .tw 


( 7 ) 


If a jet breakup is assumed according to Rayleighs theory (linear), 
the drops will be separated by equal distances 

A* = 4.51 Dj 


and 


■d = 1.89 D, 

u 
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Then the time interval between two drops will be 


,*■ 


t. - t. 
i i-i 


U 


{ 8 } 


and the ratio of the pulse interval (= time interval) and the pulse 
width 


X* / d 
U / U 


= constant 


(9) 
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Appendix 2 


Estimation of the initial disturbance amplitude 

From the linear theory of the jet breakup (Rayleigh, 1878) we find for 
the jet length 



where Lg is the breakup length, U, the jet velocity, D. the jet diameter, 
u the growth rate for the most unstable disturbance and So the initial 
amplitude of the disturbance. Most experiments with laminar jets show 

that In has a value of the order of 12. 

260 

The growth rate of the most unstable disturbance is found to be 

Wj / 

with a = 73.10 -3 N/m 
p = 1000 kg/m 3 
D, = 4.10“ 3 m 

tJ 

we find y* - 34 s- 1 

In our experiments we have 

Oj = 1,25 m/s and Lg - 0.19 m 


so that 



260 0 . 

J 


In 


5 







- 43 - 


_5 

and So - 1.10 m 

_8 

For a comparable laminar jet we find Sq = 1.10 m. Thus the initial 
amplitude of the disturbance for a turbulent jet is about 1000 times 
that for a laminar jet. 
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SYMBOLS 


c excentricity c =\/l-b 2 /a 2 , a arid b the major and minor axis 

cf calibration factor 

d drop diameter 

f frequency; f* = O./A* 

vJ 

h pulse height 

1 width of the light field interruption 

n number of samples 

t time 

tw pulse width 


D tube diameter; D., drop diameter; D-, jet diameter. 

U J 

H dimensionless pulse height 

P dimensionless pulse interval; p = jP, j = 1,2,3 _ n 

Q flowrate 

U Eulerian velocity; D, time mean value; u, turbulence component; 

u 1 =\/u^, root mean square turbulence velocity component; sub¬ 
scripts j, t refer to jet and tube respectively. 

Vj drop volume 


x 

v 

a 

P 

T 


wavelength; x* , most unstable wavelength 

kinematic viscosity 

surface tension 

density; p , ambient density 

dimensionless time 


Re 

We 

Sr 


Reynolds number; Re^ = 

Weber number based on the ambient density; We 

f 0. 

Strouhal number; Sr = -- 


p a D J°i 







typical breakup curve 


figure I. 1 
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figure II. 1 





































































A. Open tank 

B. Overflow system 

C. Direct connection to the normal 
water supply 

D. Controll valve 

E. Flowmeter (Fisher and Porter) 
y.Jx: -JV: '''-j\ F. Calming chamber 

G. Test section 


figure II. 2 
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figure II. 3 
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figure II. 4 
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fossa navicularis with elliptical meatus 



Profile Y-Y 



Profile Z-Z 


figure II. 5 

























































A. light source (Philips 7027 15V50W) 

B. 1ens ; doublet. 

C. D. slits 

E. photodiode (U.D.T. pin-L4) 

F. amplifier (first stage) 

G. drop train 


LABORATORY DROP SPECTROMETER 


figure II. 6 
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Figure II. 7 



























































L.D'.S. signal (digital) 
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figure IV. 1 
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figure IV. 3 



Restriction 5 placed at several distances from the 


tube orifice; Rep = 5000 



figure IV. 4 





Restriction 5 placed at several distances from the 
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Mean value and standard deviation 


(72 groups 


1024 
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-15 sec duration) 


case 
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Restriction 5 placed at a distance of 10, 

20, 30, 40 er 50 diameters upstream from 

the orifice 

(= case 1, etc. 5) 


Table IV. 2 
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Different rigid restrictions placed at a dis¬ 
tance of 10 diameters upstream from the tube 
orifice 

open tube 
A restriction 1 
•t restriction 2 
X restriction 3 
O restriction 4 
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figure IV. 7 
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Different restrictions placed at a distance of 10 and 40 diameters upstream 
from the tube orifice (open tube = case 0, restriction 1 = case 1, etc) 
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Different restrictions placed at a distance 
of 10 diameters upstream from the orifice 

(open tube - case 0; restriction 1 = case 1 
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Different restrictions placed at a distance 
of 40 diameters upstream from the orifice 

(open tube = case 0; restriction 1 = case l,etc) 


Table IV. 5 
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